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SUMMARY 


The application of NASTRAN to the structural analysis of the Tokamak Fu- 
sion Test Reactor (TFTR) Toroidal Field (TF) magnetic coils and their support- 
ing structures are described in this paper. The primary applied loads on the 
TF coils are electromagnetic and thermal. The complex structure and the tre- 
mendous applied loads necessitated computer type of solutions for the design 
problems. In the early stage of the TF coil design, many simplified finite 
element models were developed for the pvirpose of investigating the effects of 
material properties, supporting schemes, and coil case material on the stress 
levels in the case and in the copper coil. In the more sophisticated models 
that followed the parametric and scoping studies, the isoparametric elements, 
such as QUAD^l , HEX8 and HEXA, were used. The analysis results from using these 
finite element models and the NASTRAN System were considered accurate enough to 
provide timely design information. These analytical results were further con- 
firmed by using the EBASCO developed isoparametric composite element. 


INTRODUCTION 


The TFTR, now under construction and expected to be operational by I 98 I, 
will be the first Tokamak in the world to produce l^J MeV fusion neutrons and 
3.5 MeV fusion alpha particle? in a break-even experiment involving a deuter- 
ium-tritium plasma. The information to be obtained from the TFTR project will 
contribute to the effort aimed at the construction of demonstration fusion pow- 
er reactors in the ea.rly 1990s . This would help meet our nation ' s energy de- 
mands . 

The main components of the TFTR are the large doughnut- shaped vacuum ves- 
sel (or torus) which contains the plasma and the TF coils which generate the 
magnetic field for plasma confinement (see Figure 1). There are 20 TF coils 
(see Figure 2); each coil has hk turns of copper winding which are insulated 
from one another. The coils are wrapped in epoxy and placed in a structural 
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case (see Figure 3). These coils Eire supported near the machine's center 
(see Figure 1) hy an inner ring and column system and outboard of the coils 
by outer shear -compression members (see Figure 4). 

The coils are operated in a pulsed mode with the current increased from 
zero until the desired magnetic field strength is reached. The field is held 
constant (flat top) for approximately 3 seconds before decreasing. The mini- 
miun time between pulses is 300 seconds; the coils are required to sustain 
300,000 full -power pulses. Heat generated by the electric current in the cop- 
per coils is removed by cooling water flowing through the coils at an inlet 
temperature of 10°C (50°F) and an outlet temperature not more than 68.3°C 
(155°F). 

The primary applied loads on the TF coils are electromagnetic and ther- 
mal. There are two electromagnetic loads; one is due to the interaction of 
the toroidal current and its field which results in a net centering force of 
27 076 hK (6087 kips). The other electromagnetic load is due to the interac- 
tion of the toroidal current and the Equilibrium Field (EF) which results in 
a lateral moment of 11 890.O kNra (105 24o inch-kips) (s e Figure 5). The 
thermally- induced load is mainly due to temperature difference between the 
copper winding and the coil case. 

The complex structure and large applied loads necessitated computer-based 
analysis for solutions to various design problems. In order to investigate 
the effects of material properties, supporting schemes, coil case material and 
manufacturing tolerances on the stress levels in the case and in the coil, 
simplified finite element models were developed and analyzed using NASTRAU. 

In the more sophisticated models that followed the parametric and scoping 
studies, isoparametric elements such as QUAD4, HEX8 and HEXA were used. 

This paper will describe the various finite element models that were used 
in support of the coil design. In addition, the element types, size of each 
model, multi-point constraints and cyclic symmetry features of NASTRAW are 
discussed. The analytical results and correlations among the various models 
are also presented. A comparison of the results obtained using the common 
HASTRAN isoparametric elements with results obtained using the EBASCO-devel- 
oped composite isoparametric element will be included. 


FINITE ELEMENT MODELS 


The finite element method was employed in the structural analysis of the 
TF coils. Several models of the TF coil were generated during the preliminary 
and the detailed design phases of the Tokamak Project. They can be catego- 
rized as follows: 

• Simplified 2-D and 3-D Models 

• Detailed 3-D Model 

• Detailed Local Models 

• Detailed 3-D EBASCO developed Composite Element Model. 
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The first group of finite element models (FEM) vas used primarily during 
the preliminary design phase, vhere numerous studies were conducted to achieve 
an optimum TP coil design with respect to coil shape and coil support schemes. 
Since the prellminai'y design phase i;' character iaed by fast design changes, 
the PEMs must be si‘pl“ enough for modifications to be accomplished in a short 
time and the results generated must meet the design objectives. The first 
three PEMs, summarized in Table I, are considered as Simplified 2-D and 3-D 
models. TF1G3D was later modified to include dynamic, fault and sei.gmic anal- 
yses . 


Following the completion of the design studies using the Simplified Mod- 
els, a preliminary design configuration of the TP coil assembly was subjected 
to a more detailed stress analysis. The second generation 3-D model (TF2G3D 
in Table l) was developed for this effort. This model includes the TP coil, 
outboard inter -coil shear’ compression box and the inner rings attached to the 
central column. Analytical results from this model had shown that the stresses 
in the epoxy and the copper winding in the vicinity of the ring supports were 
high, but this model was not fine enough to provide accvrrate results required. 
Since the demand for design information in the epoxy, copper, coil case bolts 
and the bolts in the yoke were of such a detailed nature, Detailed Local Mod- 
els were developed. These models are suinraoriaed in Table I as TF3G2D, TF3G3D, 
90 ° Yoke Segment, and Local Stress Concentration. 

These models provided timely information for the TP coil design effort 
from the preliialnai'y design to a final design phase of the Tokamak Piograjn. 

In the meantime, a more sophisticated isoparametric composite element was being 
developed. A model using the composite element was generated to confirm the 
results produced by previous models using available KASTRAH elements, as well 
as to verify the structural integrity of the TP coil. 

A more detailed description of the above mentioned PEMs are presented in 
the following paragraphs. 


First Generation 2-D Model (TF1G2D) 

The purpose of this model was to study the different coil shapes and to 
inves Igate various support conditions. 

The first generation 2-D model consisted of tlrree differtnt models; the 
circular shape, the dee shape and the tear drop shape, which was basically the 
dee shape turned around. Due to the synanetric condition of the geometry, as 
well as the loading (TF force), only one-half the coil was considered in this 
model, along with its supports. 

The coil was represented by a one-line model consisting of beam elements , 
with the nose and outboard supports simulated by spring elements (see Figure 
6). All of the beam elements had six degrees of freedom per element, while 
the springs had two degrees of freedom per element. 

For all analyses, both homogeneous and composite section properties were 
considered, along with the change of the cross section due to the wedge at the 
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nose. These properties were calculated at the center of each member, while 
also applying a factor of the ratios of the Young's moduli of the steel a.nd 
copper, in order to compensate for the different materials. 

The model had support points along the horizontal line of symmetry. 

These points were assumed to only have the freedom of horizontal displacement. 
The supports at the outboard inter-coil connections were simulated as horizon- 
tal springs whose stiffnesses were calculated based upon the resist • .ce capa- 
bility of the connections in the horizontal direction. At the nose, the coil 
originally had a wedge support where the side plates of the ?teel casing actu- 
ally came in contact with the adjacent coils, providing resistance to horizon- 
tal displacement. In later design, this was replaced by support rings coming 
from a center column. All three models were considered with the wedge support, 
and in addition, the circular shape model was analyzed - ith the support rings. 
The values of these spring constants depended on the stiffness of such support. 


Second Generation 2-p Model (TF2G2D) 

The purpose of this model was to understand the force and moment trans- 
mitted from the copper winding, through the epoxy, to the steel case, as the 
coil was subjected to the toroidal field force. This study included the ef- 
fects of shear modulus in epoxy, the Young's modulus of the copper winding, and 
the failure of the epoxy in shear. 

Due to the symmetric condition of the geometry, as well as the loading 
(TF force), one-half of the coil with the supports was considered in this 
model . 

The model consisted of three radii, with the inner and outer ones contain- 
ing the steel casing grid points and the middle one for the copper grid points. 
When the constant copper cross-section design was introduced, another radius, 
coincident with the copper grid point radius, was added for additional steel 
casing grid points. This was needed in order to create a finer mesh to handle 
the varying thickness of the side plates. 

In the model, the steel casing, including cuter and inner flanges, and 
the side plates were modeled as BAR (beam) and QDMEM2 ( q.uadrilateral membrane) 
elements, respectively. The copper winding was modeled as BAR elements which 
connected the copper nodal points located in the middle of the coil. The ep- 
oxy between steel casing and the copper winding were represented by ROD (beam 
with pin connections) and SHEAR elements. The cross-sectional areas of ROD 
and SHEAR elements were calculated based on the equivalent stiffness of the 
epoxy in resisting the relative movement of the steel casing and the copper 
winding in radial and tangential directions, respectively. The supports at 
the nose and the outboard box of the coil basically resisted the movement of 
the steel casing in the horizontal direction and were modeled as horizontal 
springs (ELAS2). 

Many cases have been investigated with this model and various support 
schemes have been analyzed. The bond in the epoxy between the copper bundle 
and the casing have been considered, as well as the case of nc bond in epoxy. 
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Analysis has also been pei‘foi’iiied to determine the redistribution of stresses, 
moments, and forces with varying .uaterlal and cross-sectional properties. 

This model has also been used to determine the effect of having a gap netween 
the coil and support rings. In addition, various analyses have been performed 
using the constant copper cross-section design. 

The results obtained from these analyses included displacements of grid 
points, loads in the copper, forces at the supports, and the forces and 
stresses in the steel casing, copper, and epoxy. All of these results were 
obtained for both toroidal field (TF) loads and for thermal loads. 


First Generation 3-D Model (TF1G3D) 

This model, shown in Figure 7, was originally generated to investigate 
support schemes for the TF coil subjected to out-of-plane loading. It was 
later modified for dynamic (EF out-of-plane), favilted (TF coil short circuits) 
and seismic analyses. This coarse type model was selected to avoid the expen- 
sive computer cost usually incurred in analyses for dynamic and nonsynmietric 
static loads. It was estimated that the effects of these loads are minor on 
the stress in the criticai region of the TF coils and the inner support struc- 
ture (ISG) compared with the normal operating electromagnetic static loads. 

This model can give a qualitative comparison of TF coil and ISS stress results, 
and reasonably accui’ate results for reaction forces at the supxjort pedestals 
and floor system. This section first describes the finite element model used 
in the nonlinear dynamic analysis performed for the out-of-plane TP coil pulse 
loading induced by the interaction of the EF coil magnetic field witb the TP 
coil electrical cm-rent field- and then describes the modifications made to 
this model in the seismic and faulted load analyses. 

Dynamic Analysis Model - This model has taken into account the geometric non- 
linearity associated with the TF coil pedestal supports, which prevent hori- 
zontal lateral displacement of the TP coil in only one direction. Structural 
damping was also included. 

Because of the geometric rotational syimnetry of the TF coils, ISS and 
floor system, a tlu-ee-dimensional finite element model of a 36° segment of the 
structure, composed of one-dimensional and two-dimensional finite elements, 
was used in the analysis (see Figiu-e 7)- The model was sufficiently detailed 
so that the effects of the TP coil pulse loading, seismic loads and fault 
loads coiild be determined, but local structural details, such as keyed and 
bolted splices and cutouts, were not incorporated into the model. A total of 
262 nodal poincs and 248 finite elements (114qUAD4, 76 BAR, 44rBAK, 12E0D, 
2CELAS2) were used. 

The TP coils, represented by BAR elements, were connected to the support 
pedestals hy means of REAR elements and CELAS2 spring elements which talce into 
account the coil base pedestal stiffness. These connections provide horizon- 
tal lateral load transfer from the TF coils to the pedestals. At these con- 
nection points between coil base and pedestal top, a nonlinear tx'ansient 
force-displacement relationship was defined to represent the nonlinear behav- 
ior of the pedestal support system. This geometric nonlinearity is such that 
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the restraint imposed on the coil base is a function of the direction of the 
relative horiaontal lateral displacement between the coil hase and the pedes- 
tal. Physically, this restraint is provided by a snubber plate against which 
one side of the coil base initially rests. The snubber plate location alter- 
nates from side to side, from coil to coil. Under load, when the coil base 
tends to displace such that it bears against the snubber, resistance to that 
displacement is provided. When the coil base displaces away from the snubber, 
it can do so freely and without restraint. 

The floor structure was modelled as a system of QUAD'+ elements and 
elements. Connection of the pedestals to the midheight of the floor were made 
w'th RMR elements. 

The TFTR structure is externally supported at the machine floor level by 
connection of the outer edge of the machine floor to vertical columns and the 
building floor system, and by connection of the inner edge of the floor system 
to vertical coliunns. These external boundary conditions were incorporated 
into the finite element model by fully restraining all degrees of freedom at 
the outer edge of the machine floor (3 displacements and 3 rotations). The 
vertical floor columns were not incorporated into the model because their ef- 
fect on the analyses performed is negligible. 

To simulate the symmetric behavior of the structure, boundary condition 
relationships equating displacements and rotations at boiuidary grid points 
were imposed. 

Model Modifications for Seismic and Faulted Load Analyses - The primary modi- 
fications made to the model for the seismic and faul^ ’’ \o, d analyses were 
related to the boundary con 'itions. Since these lo^ were not rotationally 

syimnetric like the dynamic load, the syiimietric bounds . condition used in the 
dynamic analysis was not valid. KASTRAN's cyclic symmetry featui’es were used 
to represent the entire 360° structure, including 20 TF coils. In order to 
use the cyclic symmetric features, the non-linear connections between the TF 
coils and the tops of the pedestals were removed. It was assumed that only 
one out of every pair of coils was connected to the support pedestal . The 
floor system, which has a minor effect on these analyses, was totally re- 
strained to reduce the problem size. 


Second Generation 3-D Model (TF2G3D) 

The first generation finite element models described in previous para- 
graphs were used mainly in the preliminary design of the TF coil, in the selec- 
tion of coil shapes, and finally in the baseline design. Having designed a 
baseline TF coil, the next step was to perform more detailed stress analyses. 

In order to support the detailed design phase of the TP coil, a more sophisti- 
cated finite element model wan required. The second generation 3-D model was 
developed for this effort. This model is a thi’ee-diinensional finite element 
model which includes the TF coil, outboard inter-coil support box and the 
i.i»jv.r rings attached to the central column, support. The second generation 3-D 
finite element model represents an l8° sector which includes part cr the ring 
and column support, one coil, and two half sheai' -compress ion supports from 


252 



eAoli aide ol' the coil {see Figui’e 3). The coil case is simulated by plate 
elements (QUAD2, TRIA2) with inplane and bending stiffnesses. 

There are J|!| copper windings in each coil; they are electrically isolated 
from eacli other and the coil case. The copper windings are represented by 
four solid elements (HEX8), and they ai'e sepai*ated from the coil case by the 
epoxy and the electrical insulating material. 

The center column and ring suf'po'.'ts and the outbooa'd inter-coil sliear- 
compression boxes are represented by plate elements. The TF coil finite ele- 
ment model is supported at the bottom of the coil. Tills pedestal support al- 
lows the coil to move in the upward direction (no hold-dowi) but supports the 
weight and loads in the dowiwai'd direction. In addition, the coil is oJ.lowed 
to move laterally (pex’pendiculai’ to the coil plane), but is restrained by a 
key in the opposite direction. Radially from the device center, the coil is 
unrestrained, The coil nose oi'ea is connected in three directions to a 4- 
rihg system. The two outer rings ai-e in tiu’ii attached to the central cylin- 
der only in the tangential direction, nox'iiial to the coil plane. The inner 
rings, however, are rigidly connected to the central cylinders. 

The cyclic syimiietry capability of NASTRAW was used to effect r..;e remain- 
ing 19 sections of the TF coil system. The boundai'j'' points ai'e defined as in- 
j5ut to the program. The computer time was excessive when cyclic si'iiimetry was 
used. Since the sti'uctiu'al analysis of the baseline configuration was mainly 
using sjnnmetric load, the multipoint constraints (MFC) can be used on the 
ooundary points and, therefore, one coil model is needed. The computer time 
has been reduced considerably by using the Nff*Cs. However, for unsyimuetric 
loading, the cyclic symmetry must be used. 

The epoxy bond between the copper and the coil case may have been sepa- 
rated during the curing process. In oi'der to simxilate this condition in the 
finite element model, the epoxy can have only compressive strength. Wxen the 
resxilt shewed any epoxy spring element that was in tension, the spring element 
was eliminated fi-om the model. After several iterations, all spring elements 
that were in tension were eliminated from the model. The results were repre- 
sentative of the no epoxy bond case. 

The second generation 3-D TF coil finite element model has IOO 6 grid 
points, 127 I 1 elements (1|6 TR1A2, 524 QUADS, 1T6 HEX8, 523 ELAS9), and 4239 de- 
grees of freedom. A typical run of this model, employing MFC boundary condi- 
tions, had vised approximately 3000 system resource vmits (SRU) of computer 
time, whereas a run using cyclic syiimietry had used three times as much time. 

There v/ere several TJASTRAlf rvuis made using the 3-2 t'lodel in direct sup- 
port of the TF coil design effort. The ex'fect of different materials (e.g., 
titanium oi‘ steel) for the coil case and the inner support strvictvu’e was in- 
vestigated. Stvidies were done to vary the location and nvunber of rings for 
the central colvuiin support. Vai'yiag thicknesses of the coil outer flange was 
also analyzed using the 3-D model. 

The grid point forces and displacements, element forces, and stresses are 
pi'inted output from the NASTRAE run. Some of the results ore usually svmmia- 
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I'lzed manually to provide immediate design information. However, the amount 
of data remaining to he siumnarized was massive; therefore, it was imperative 
that a computer aided post-processing program he developed. The output of 
these FORTRAN programs was in both tabulated and plotted formats. Plots of 
shear load per unit length for the inner and outer flanges were made. These 
plots provided the designer with the necessary infoinnation for sizing the holts 
required for fastening the coil sideplate to the inner and outer flanges. 

These plots and tabulated data were subsequently used for fatigue analysis. In 
addition to these plots, copper hoop stress, epoxy normal stress, and casing 
membrane stress intensity were also plotted. 

The copper stresses were printed in the basic rectangulai’ coordinate sys- 
tem; however, this was not convenient for evaluation. The latest version of 
KASTRAN with the HEXA element has rectified this shortcoming. However, the 
corner stress "Jiuup" for the solid element was very disturbing. Fortunately, 
the copper hoop stress was reasonably acciu'ate when the Ij-corner HEXA stresses 
were arithmetically averaged. This was perhaps due to the fact that there was 
no applied load in that direction. In the radial direction, the HEXA stresses 
cannot be simply "averaged". It was estimated using the grid point force out- 
put. More experience is required for the interpretation of the solid element 
stress output. 


Third Generation 2-D Model (TF3G2D) 

Since the inter-laminar epoxy stresses from the second generation 3-D 
model were obtained from a relatively crude epoxy/copper simulation, a more 
detailed third generation 2-D model (TF3G2D in Table l) was developed. Sym- 
metry required only half of the coil in the model (see Figure 9)- Bar ele- 
ments were used for copper layers and inner and outer flanges of the case. A 
shear panel provided only shearing strength for the copper -to-copper and cop- 
per-to-case interfaces. A rod element provided two degrees of freedom which 
gave direct relative motion of the interfaces. An energy equivalence was made 
to obtain a spring constant for the rod and a shear strength for the shear 
panel. In addition to the elements in the bar-shear-spring model, a membrane 
element iras applied in TF3G2D. This type of element was used to simulate the 
in-plane deformation and strength of the case side plates. 

Local bending and tangential extension of all grid points at both ends of 
the half coil were restrained due to symmetry. Rod elements provided horizon- 
tal elastic supports at center column ring suppoi-ts and connecting box sup- 
ports, An energy equivalent was made for those rod elements to obtain the 
proper stiffness constants. 

The number of nodal points for a representative model of this kind was 
612, which resulted in l6o^t degrees of freedom. The two-dimensional elements 
included 900 bars, 480 rods, 120 membranes, and 390 shear panels. Total niam- 
ber of force units applied to the system, for magnetic force in static analy- 
sis alone, was 434, Number of variable temperature grids was 595. Average 
computing time for each rim was 648 system seconds on a CDC-6600 for tliree 
loading cases. 
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Third Generation 3-D Local Model (TF3G3D) 

The third generation 2-D model described in the previous paragraphs was 
used to investigate shear stress distribution in epoxy due to in-plane loads 
(TF and thermal). The third generation 3-D local model (see Figui'e 10) was 
used primarily for the study of the sheeir stress in the epoxy and copper 
stress due to in-plane and out-of-plane (EF) loads at the critical area. 

To keep the model small enough to be analyzed in the current computer and 
still render good critical stress data, a portion of the TF coil was analyzed. 
This local model consisted of the coil cut from horizontal zero degree plane 
and approximately 5^i-degree plane from the horizontal. By applying a symmet- 
ric and an anti-syimiietric boundary condition to the model, a cut was made to 
the mlddJ.e plane which was the separation of the reverse copper winding by a 
barrier strip. The barrier strip retained full capability of bonding; how- 
ever, the bonding between the copper bundle and the case was removed except 
where bhe compression A'as developed by the combined loads. The nose support- 
ing ring was simulated by springs in all directions. Bi" undary forces were ap- 
plied to the 5^1-degree cross section. These forces were obtained thi’ough 
global 3-D model. Magnetic loads and thermal data were carefully distributed 
to all grid points oi the model. It was the thinner copper layer and the 
flexibility of the layer in the other direction that made the usage of plate 
elements in 3-D model. TF load applied to the model included a radial in- 
plane load and a transverse out-of-plane load which sq.ueezed the copper to the 
barrier strip. 

The differences between this model and the third generation 2-D model 
were that plate elements were used in the 3-D model to represent copper lay- 
ers, whereas bar elemerts were used in the 2-D model. In addition, solid el- 
ements were used in the 3-D model for the outer flange, whereas bar elements 
were used in the 2-D model. Due to the geometric complexity and the thickness 
variation of the outer flange, the solid element was used. Two types of solid 
elements were used in the 3-D model. A HEXA element had eight nodal points at 
the corners. Each corner had three translational degrees of freedom. Two 
HEXA elements used in the model had extra points along the edges for the 
proper modeling of wedge transition. A WEDGE element had six nodal points at 
the corners. Each corner had tlu‘ee translational degrees of freedom. HEXA. 
elements were used also for the boi-rier strip. Plate elements used for cop- 
per, side wall, and inner flange were QUADUs. A QUADU element had four nodal 
points at the corners and each nodal point had five degrees of freedom. The 
cooling hole stiffness was reduced by using offset capability of QUADh . Shear 
panels used in the model were SHEAR'S which were obtained from energy equiva- 
lent. The shear panels had only shear capability and each panel had four de- 
grees of freedom. These panels connected the copper layers in two ways. A 
series of shear panels took the resistance of slip between copper layers in 
tangential direction and the other series of shear panels prevented copper 
from relative movement in the transverse out-of-plane direction. Finally, a 
set of two degrees of freedom rod e.-oments provided compressibility of copper 
layers and epoxy which wrapped copper layers and copper bundle. The connec- 
tion between the outer flange and the side wall was accomplished with bars 
and rods. 
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The nuraher of grid points in the system was 2,3^7- had 10,T61 de- 
grees of freedom for symmetric load condition and 11,231 degrees of freedom 
for anti-symmetric load condition. The number of total bulk data lines (or 
cards) was lH,89^. Total execution time in SRU was lh,6l0 on a CDC- 66 OO for 
three subcases and two subcombs by NASTRM code. 


TF Coil Yoke Finite Element Model 

The yoke finite element model is a 90° segment (90° - l80°) taken from 
the second generation 3-D model. The boundary loads were based on the TP2G3D 
model. The main objective of the 90° segment yoke model was to ascertain the 
bolt loads between the yoke and the coil case, and between the yoke and the 
shear-compression box (see Figure 11). Revisions to the yoke model required 
fewer manhours and computer usage time compared to the full coil model. 
Therefore, this model was used for yoke design support and parametric studies. 
The yoke model includes all elements associated with the second quadrant 
( 90 ° - 180 °) of the second generation 3-D TF coil model. The coil case is 
simulated by plate elements (QUADt, TRIA.3) with inplane and bending stiff- 
nesses. The copper windings are represented by 6 solid elements (HEXA) ; they 
are separated from the coil case by the epoxy and the electrical insulating 
material. Due to design considerations, the epoxy material between the cop- 
per winding and the casing is assumed to take compression loads only, and has 
no shear or tension capabilities. The spring element (ELAS2), which has two 
degrees of freedom, is used to simulate these compression forces. 

The bolt shear and contact stiffnesses are simulated by using the spring 
element (ELAS2). The yoke and shear-compression box are simulated by plate 
elements. A typical section is shown in Figure 11. The bolt shear stiff- 
nesses were calculated by assuming that the parent material had yielded to 1 % 
of the bolt diameter, and these stiffnesses were beamed to the nearest grid 
points . 

The coil case is made of steel; the gages at the side wall, outer and 
inner flanges are 2.9 cm (l 1/8 in. ), lit .6 cm ( 5 3/i^ in. ) , 6 .9 cm ( 2 3/it in. ) , re- 
spectively. The shear-compression box and yoke assembly are shown in Figure 11. 

The yoke model uses the WASTEAE multipoint constraints (MPC) at the 
shear-compression box boundary points. The boundary loads at both 90° and 
180 ° are calculated using the TF2G3D model. 

The contact stiffness between the shear -compression box and nhe yoke and 
the bolt axial stiffness are simulated by using spring elements. When the 
contact spring elements are in tension, they are eliminated from the input 
data set. When the bolt is in tension, the tensile stiffness is used. 

The yoke finite element model has 616 grid points, 922 elements ana 
3,696 degrees of freedom. The computer time usage is 1,500 seconds (SRU). 
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Local Sti-ess Concentration Model 


The fatigue life requirement is one of the critical items in the TP coil 
design. A better understanding of the peak stress (including stress concen- 
tration) distribution around the bolt hole of the coil casing is essential to 
the estimation of the fatigue life. The primal” contributors to the peak 
stress aroimd the bolt hole are the bypassing stress (hoop tension) and the 
pin loads (both tangential and radial directions). Because of the complicated 
natiu'e of the casing geomebry, the effects of the by-passing stress and the 
pin loads are coupled. It implies that without knowing the load path around 
the bolt hole, it is difficult to predict how much of the peak stress is con- 
' ■ ibuted by the by-passing stress or by the pin loads. Therefore, a finite 
(.'leraent analysis was initiated to study the load path in the critical bolt 
hole and predict the critical peak stresa. The finite element model included 
a segment of side plate with the critical bolt holes connecting the inner 
flange , 

The "Multilevel Superelements" approach available in the NASTRAN progx'am 
was cl; 'Sen in this analysis. The advantage of the superelement approach is 
to allow a fine mesh model near the critical bolt holes, but a coarse mesh 
away fi’om the critical area without losing the accm-acy. Moreover, it is only 
required to have one fine mesh bolt hole model which can be repeated as many 
times as desired to simulate a series of bolt holes. 

Two mesh generators were written to automatically generate two separate 
models. One is the coai’se mesh model, designated as the residual superelement 
(see Figure 12), which included a segment of the side plate, including the 
bolt holes away from the critical area. Tlxe other is the fine mesh model of a 
typical bolt hole, designated as the px’irnary superelement (see Figure 13). 

Then any number of the secondary superelement, which is similar to the primary 
superelement, can be automatically generated. Figm'e lU shows the combined 
geometry of the residual and primai’y superelement s , and indicates the location 
of two secondax’y supe 'elements . The mesh generators were programed to include 
the parameters, such as diameter of bolt hole, pitch of bolt hole, numbers of 
bolt holes and the dimensions of the side plate segment . 

The model is two-dimensional. The quadx’ilateral (QUADIi) and triangular 
(TRMEN) plate elements wei*e used in the residual superelement, while the 
TRMEN element was used in the pi'imai’y superelement . The QUADl-! element has 
eight degrees of freedom, while TRMEN has six degrees of freedom. 

The force boundary condition was used to describe the boundary at the 
edges of the side plate segment and the pin loads acting on the bolt bole. 
These data were obtained from the resuJ,t of 3-L analysis (TF2G3D model). The 
residual superelement model was vised to check these applied forces for the 
static eqviilibrium. The stresses in the side plate resulting from the equi- 
librium check run were also compared xn-th the 3-D result. The unifox’m radial 
pressure distribution was assumed for the pin loads acting at the px'imax'y and 
secondary superelements. 

Three different loading conditions were considex'ed fox' each model; pin 
loads in both tangential and radial conditions, pin load in tangential direc- 
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tion only (sheeir key to take radial shear), and no pin load. The "no pin 
load" condition was not the real design option being considered. However, it 
included to approximately simulate the condition if the friction took all 
the shear loads. For all loading cases, the stress printouts plus the stress 
contour plots for the maximum principal stress were obtained directly from 
NASTRAH runs . 


Composite ElemenL 3-D Model 

This model uses the EBASCO developed 3-D solid composite elements for 
heterogeneous material. These elements simulate the composite action of the 
copper and epoxy bundle and also give detail stress information in both cop- 
per and epoxy. The composite element is defined by 8 to 20 nodes. 

The coil case was represented by composite elements with homogeneous 
isotropic properties. There were two radial layers to represent the 5 3A 
inch thick outer flange. The yoke, inter-coil shear compression boxes, pedes- 
tal support, and support rings at the nose were also modeled using the com- 
posite element. The composite element TF coil model has 25^7 grid points, 

T^tO composite elements, and 736 non-linear spring elements. Further discus- 
sion and theoretical presentation of the composite element can be obtained 
from references 1 and 2, respectively. 


DISCUSSION OF MODEL RESULTS 


Representative results from the previously described finite element mod- 
els are discussed in the following paragraphs. 

Figure 6 shows the various coil shapes represented by the TF1G2D models. 
The high energy requirements for the dee shape had eliminated it from further 
consideration. The practical design diffic^LLties of the inclined spring sup- 
port for the tear-drop shape with T1 = T2 (see Figure 6), and the large ten- 
sion jump which exists in the tear-drop shape with T1 ^ T2 have also elimin- 
ated it from consideration. Manufacturing simplicity of the circular coil and 
its low-energy requirements were the important points for the selection of the 
circular configuration. 

Figure 15 shows the casing stress versus location and stiffness of the 
double ring supports for the circular coil. This was obtained by using the 
TF2G2D model. 

Using the TP1G3D model, a NASTRAN non-linear dynamic transient analysis 
was performed for the 0.42 Tesla Strong Compression plasma off run 4l EF load 
condition. This out-of-plane TF coil loading was applied to TF coil nodes as 
concentrated time dependent nodal forces. The load-time profile used is' 
shown in the normalized strong compression EP current of Figure l6. Both 
dynamic and static analyses were made using the TF1G3D model. Comparisons of 
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the ytressea obtained from tlie dynamic and static analyses were made to deter- 
mine the dyncauic effect of the operational EF loading. Typical stress multi- 
plication factors for various parts of the TF coil structure were less than 
1,09 except in the area near the pedestal support. The factor is l.l8 at the 
pedestal support. The tangential load transmitted to the floor at the pedes- 
tal increases from 0 (static loads) to 60 kips (dynmuic load). This is due 
to the Inertia load induced by the restraint at the pedestal. Similar’ re- 
sults were obtained for the faulted load (one or more coil short-circuited) 
and seismic analyses. Both the faulted and seismic loads have little effect 
on the coil, whereas the pedestal receives a significant amount of loads. 

A plot of maximum shear stress in the epoxy due to the TF and thermal 
loads from the TF3G2D model is sh.oim in Figui’e IT. Since the model is 2-D, 
epoxy shear stress due to EF load is not available; it can be obtained from 
the 3-D model TF2G3D (see Figvu’e l8) . 

The design criteria used for the coil case is basically derived from the 
ASME Boiler and Pressure Vessel Code. This code uses stress intensity, which 
is defined as the largest algebraic difference between any two of the thi'ee 
principal stresses. In order to facilitate the evaluation of coil case stress, 
all NASTRAN output for the TF2G3D model were post-processed and presented in 
both tabulated and plotted formats. A typical plot of membrane stress inten- 
sity is shoTO in Figure 19. A plot shown in Figure 20 was used by the designer 
in the determination of bolt requirement (the bolts i’asten the coil side plate 
to the flange). 

Stress distributions in various parts of the TF3G3D local model are shorn 
in Figure 21. These results are for the early TF coil, with titaniiuu as cas- 
ing material. This model was not re-run for the latest configiu'ation , since 
, at that time the composite element was available for model development. How- 
ever, the results from the eai'ly coil configiu'ation served as a benclimark for 
tlie TF2G3D model, 

Tlie yoke load distribution calculated by using the yoke model is shown 
in Figure 22. 

The major principal stress and the maximum stress intensity (twice of 
inaximiuti sliear) ai'e presented in Table II. In general, the maximum principal 
stress is the tensile stress in the tangential direction of the TF coil. The 
stress level was mostly due to the by-passing stress; the radial pin load con- 
tributed about 15/a to l8^ and the tangential pin load contributed about 7 %. 

The difference between the principal stress and the stress intensity was pri- 
marily due to the bearing (compressive) stress of tlie pin loads. In load 
cases 2 and 3, there were minor or no differences between the maximum princi- 
pal stress and stress intensity, since no I'adial pin loads were directly ap- 
plied at the edge of the bolt holes. 

The stress levels in the coil components from the Composite Element mod- 
el and the second generation 3-D model are compared in Table III, Despite 
the number of different model techniques used in these models, the stress 
levels were fairly close. 
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CONCLUP^Aia REMARKS 


There were sufficient tJASTRAH d .''.(gnostics used in all models described 
in this paper such that good numerical results were obtained. Force bali.nce 
at the inner support structure and selected coil nodal force balances were 
perforfned so that the confidence in these models was ascertained. Further- 
more, the stress outputs from these models were compared and were found to be 
relatively close. The analysis results from these finite element models and 
the NASTRAK System were considered accurate enough to provide timely design 
information. 

The application of available MSTRAN elements to the ToFamak toroidal 
field coil finite element models was considered very effective and the anal- 
ysis results were acceptable. However, care should be taken in using the 
spring and isoparametric elements. Some of the comments pertaining to these 
elements are isted below: 

• The plate element has no inplane rotational stiffness. Special care 
should be taken when it is connected to one-dimensional elements. 

• The spring element used for non-colinear connection will induce un- 
balance moment. 

• The nodal force balance output does not include the "lock-in" force 
induced by the thermal load. 

• In a coarse model, the corner stresses for solid elements HEX8 and 
HEXA are not reliable. 

The effect of the local detail design, such as the bolt joint connection 
between the sideplate and the flange, was neglected in the TF2G3D and other 
simplified models described in this paper. This effect will be analyzed sep- 
arately when design detail is available. 
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EBASCO LOCAL STRESS q^vnifP 

COMPOSITE CQNCEN- ccpmcmt TF3G3D TF3G2D TF2G3D TF1G3D TF2G2D TF1G2D 

ELEMENT TRATION SEGMENT 


MODEL 




Table I Finite Element Models 


OBJECTIVES 

MODEL DESCRIPTION 

LOADS 

• COILSHAPESELECTION 

• COIL SUPPORT SCHEME 

• BAR, ELAS2 ELEMENTS 

• ONE-HALF OF COIL WITH SYMMETRIC 
BOUNDARY CONDITION 

TF 

• COIL SUPPORT SCHEME 
« EFFECT OF 
INSULATION 

• BAR, QDMEM2, ROD, SHEAR, ELAS2 
ELEMENTS 

• ONE-HALF OF COIL WITH SYMMETRIC 
BOUNDARY CONDITION 

TF 

• COILSUPPORTSCHEME 

• DYNAMIC, SEISMIC, 
FAULT LOADS 
ANALYSES 

• BAR,QUAD2.(QUAD4, ROD, ELAS2) 
ELEMENTS 

• TWO COILS WITH ROTATIONAL SYM- 
METRIC BOUNDARY CONDITION OR 
CYCLIC SYMMETRIC FEATURE 

TF, EF, DY- 
NAMIC, SEIS- 
MIC, FAULT 

« 

• LOAD AND STRESS IN 
CASING 

• UPDATE DESIGNS 

• PROVIDE BOUNDARY 
;C0ND FOR LOCAL 
'MODELS 

• , ESTIMATE STRESS IN- 
SIDE COIL 

• QUAD2,HEX8,TRIA2,ELAS2 ELEMENTS 

• ONE COIL WITH ROTATIONAL SYM- 
METRIC BOUNDARY CONDITION OR 
CYCLIC SYMMETRIC FEATURE 

TF,EF, 

THERMAL, 

DEAD 

• STRESS IN THE COPPER 
AND EPOXY 

• UPDATE DESIGNS 

• BAR, SHEAR, ROD, QDMEM2 ELEMENTS 

• ONE-HALF OF COIL WITH SYMMETRIC 
BOUNDARY CONDITION 

• INCLUDE TWENTY-TWO LAYERS OF 
COPPER AND EPOXY INSULATION 

TF, THERMAL 

• DETAILED'STRESSIN 
THE CASING, COPPER, 
AND INSULATION AT 
THE CRITICAL NOSE 
AREA 

• QUAD4,HEXA, WEDGE, SHEAR, ROD, BAR 
ELEMENTS 

• ONE-HALF OF 54.7“ SEGMENT WITH SYM- 
METRIC AND ANTI-SYMMETRIC BOUND- 
ARY CONDITIONS AND BOUNDARY 
FORCE FROMTF2G3D MODEL 

TF, EF, 
THERMAL 

• BOLT LOADSATTHE 
OUTBOARD SUPPORT 
(YOKE) 

• QUAD4,TRIA3,HEXA, ELAS2 ELEMENTS 

• 90“ SEGMENTWITH THE BOUNDARY 
FORCE FRQMTF2G3D MODEL 

TF, EF, 
THERMAL 

• PEAK STRESS AT THE 
BOLTHOLEOFSIDE 
WALL 

• QUAD4,TRMEN ELEMENTS 

• APPROX 30“ SEGMENT OF SIDE WALL 
WITH bom HOLES MODELED 

• MULTIPLE LEVEL SUPERELEMENT 
APPROACH USED 

TF, EF, 
THERMAL 

• CONFIRMS PREVIOUS 
MODEL RESULTS 

• VERIFICATION OF 
COIL STRUCTURAL 
INTEGRITY 

• EBASCO DEVELOPED COMPOSITE, SPRING 
AND PLATE ELEMENTS 

TF, EF, 
THERMAL 
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Table n Pealc Stress Near Bolt Hole 


STRESS IN MPo (KSI) 


LOAD- 

B IN. OUTER FLANGE 

SYa IN. OUTER FLANGE 

MAXIMUM 

PRINCIPAL 

STRESS 

MAXIMUM 

STRESS 

INTENSITY 

MAXIMUM 

PRINCIPAL 

STRESS 

MAXIMUM 

STRESS 

INTENSITY 

LOAD CASE 1 

1036.3 

1170.0 

706.0 

917.0 

(PIN LOAD IN TWO 
DIRECTION) 

(150.3) 

(169.7) 

(102.4) 

(133.0) 

LOAD CASE 2 

872.2 

872.2 

578.5 

637.1 

(PIN LOAD IN TANG 
DIRECTION) 

(126.5) 

(126.5) 

(83.9) 

(92.4) 

LOAD CASES 

807.4 

807.4 

526,8 

526.8 

(NO PIN LOAD) 

(117.1) 

(117.1) 

(76.4) 

(76.4) 


NOTE: BOUNDARY FORCES WERE DERIVED FOR 3-D (TF2G3D MODEL) RUNS 2C26 AND 2G33B 
FORTF + EF + THERMAL LOAD. 


Table III TFTR Coil Comparison of Maximum Stresses (Toroidal 
Field, Equilibrium Field and Thermal Loads) 


STRESS IN MPa (KSl) 


MODEL NAME 

COMPOSITE ELEMENT 
MODEL RUN CA-12A-2 

TF2Q3D MOOEL 
RUN2C33B 



11 ALL SOLID COMPOSITE 

I) PLATE ELEMENT FOR 



ELEMENTS 

CASING, SOLID ELE- 




MENTS FOR CORE 



2) COIL CONFIGURATION 

2) COIL CONFIGURATION 



12A 

12A 

MODEL CHARACTERISTICS 

3) FINE GRID MOOEL 

3) FINE GRID MODEL 



4) 4-RING SOLID ELEMENT 

4) 4-RING PLATE ELEMENT 



SUPPORTS 

SUPPORTS 



51 UNBONDED CASING AND 

5) UNBONDED CASING AND 



CORE (NON-LINEAR 

CaRE(LINEARSPRING 



SPRINti) 

WITH ITERATION) 


INNER 

252.3 

250.3 



(36.EJ 

(333) 

CASING STRESS 

OUTER 

292.3 

2B9.E 

INTENSITY 


(42.41 

(42.0) 

{MEMBRANE a 

SIDE 

2B3.4 

262.7 

BENDING) 


(41.11 

(3B.I) 

COPPER HOOP STRESS 

-137.S 

-101.4 



(-20.0) 

(-14.7) 

EPOXY SHEAR STRESS (Rd 

20.0 

23.4 



(2.3) 

(3.4) 




FROMTFZG30 




ANDTF3G2D 
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VACUUM 









ORIGINAL PAGE IS 
OF POOR QUALITY 

YOKE 


DIELECTRIC SPACER 




Figure 4 Shear-Compression Box and Yoke 
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27 70S kN 

INTERACTION OF TOROIDAL CURRENT (6229 KIP ) 
& TOROIDAL FIELD, 5.2 TESLA ^ 



♦ 

(6229 KIP) 
(A) TF LOAD 



INTERACTION OF TOROIDAL 
CURRENT & EF FIELD P-20 EF 
CONFIGURATION AT 28 KA 
TF = S.2T EF = 0.42T 



A-A 


(B) EF LOAD 


Figure 5 TF and EF Loads 
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(A) TF COIL 3-D FtNITE ELEMENT MODEL 


STEEL CASING COPPER WINDING & EPOXY 

STEEL CASE 
PLATE ELEMENT 



(B) COILSEGMENT OF 3-D MODEL 


Figure 8 Second Generation 3-D Model (TF2G3D) 
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(B) QDMEM2 AIMD COIMROD ELEMENTS AT SUPPORTS 


Figure 9 Third Generation 2-D Model (TF3G2D) 







ORIGINAL PAGE IS 
OP POOR QUALIT'/] 




Figure 13 Primary Superelement 



Figure 14 Remainder of Structure Considered SE 0 or Residual 
Structure 
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''INSIDE 

^OUTSIDE + f^lNSIDE 


Figure 15 Maximum Stress in Steel Casing vs Stiffness and Location 
of Support (Represented by Force Ratio) 




NORMALIZED FIELD CURRENT (EF) 


ORIGINAL PAGE IS 
OF POOR QUALITY 


(A) TOROIDAL FIELD COIL FORCE - TIME SCENARIO FOR INPLANE LOADS 



TIME (SEC) 

(B) TOROIDAL FIELD COIL FORCE - TIME SCENARIO FOR OUT-OF-PLANE LOADS 



Figure 16 TF and EF Scenarios 
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STRESS (MPa) 




MEMBRANE STRESS INTENSITY (KSi) MEMBRANE STRESS INTENSITY (KSI) 


ORIGINAL PAGE IS 
OP POOR QUALITY 


° i TF COIL MODEL 2C40B 

ELEMENT 2 SIDE PLATE 



TF COIL MODEL 2C40B 
ELEMENT 8 - FLANGE PLATE 
I nAn _ TF + PP + TMFRMAI 



THETA ' 

Figure 19 Typcial Membrane Stress Intensity 
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100 150 200 300 ICO 150 200 

MPa MPa 






MPa KSI 



'{A} MAXIMUM INTERLAMINAR 
SHEAR STRESS 


MPa KSI 



(B) COPPER STRESS IN 
HOOP DIRECTION 


CO 

MPa 


KSi 


200 


100 


r ^ 30 
00 

^ 20h 

< , 
cc 10 
m 


UJ 


0 10“ 20“ 30“ 40“ 50“ 
LOCATION IN ANGLE 


(C) INNER FLANGE STRESS 
IN HOOP DIRECTION 


Figure 21 Stress Plots from Third Generation 3-D Local Model 
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a B 0 i.TS/;'R 3 nB boundary 
ITF + er > mERMAL IN kN (KIPS) 

T BOLTS (0) 

S/C BOX TO YOKE 


13'IQ 

( 208 . 95 ) 


&1G9 
(1443.70) 

-Y S/C BOX BOUNDARY 
REACTION 


APPLIED LOADS 
0.4 


1017 

(431.E) 




(730.47) , 333 , 
3249 


i .(22).3) 
. 1G.S) > - 

i, f ' ' 

■7.41 »■ 

1 -*0 (2BB.9 

^ I 2" F 


M23.1) 

1" BOLTS (5) 


1.9) 

2" PINS (01 


|\ J COIL X (910,6) r' BOLTS (20) 
(152.6) 

(066. 1 1 

I" BOLTS (10/HOW 2 ROWS) 
YOKE TO OUTER FLANGE 

5.0D (2") PINS ( 8 I 

YOKE TOCOILSIDEPLATE 

2.54 CM [I") BOLTS (5) 

YOKE TO INNER FLANGE 


(257.B) 1" BOLTS (3) 

■ YOKE TO YOKE AT INNER FLANGE 


REACTION AT 1BQ‘ BOUNDARY 


UKIGINAL page is 

OE P<X)R QUALTTYi 


Figure 22 Load Distribution Yoke Finite Element Model 
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